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Abstract 

We explore the interplay between lines in the gamma-ray spectrum and LHC 
searches involving missing energy and photons. As an example, we consider 
a singlet Dirac fermion dark matter with the mediator for Fermi gamma- 
ray line at 130 GeV. A new chiral or local U{1) symmetry makes weak- 
scale dark matter natural and provides the axion or Z' gauge boson as the 
mediator connecting between dark matter and electroweak gauge bosons. In 
these models, the mediator particle can be produced in association with a 
monophoton at colliders and it produces large missing energy through the 
decays into a DM pair or ZZ, Z7 with at least one Z decaying into a neutrino 
pair. We adopt the monophoton searches with large missing energy at the 
LHC and impose the bounds on the coupling and mass of the mediator field 
in the models. We show that the parameter space of the Z' mediation model 
is already strongly constrained by the LHC 8 TeV data, whereas axion-like 
mediator bounds are relatively weak. We foresee the monophoton bounds 
on the Z' and axion mediation models at the LHC 14 TeV. 



1 Introduction 



Recently, from the Fermi LAT data [T], there has been an interesting observation of the 
gamma-ray hne at 130 GeV [2|[3], which may call for an explanation from physics beyond 
the Standard Model (SM) Although dark matter(DM) has not been found directly by 
underground experiments such as XENON i6|, it is an interesting possibility that the Fermi 
gamma-ray line is produced by the DM annihilation into monochromatic photon(s), such as a 
photon pair or Z7 or /17. Then, the DM annihilation cross section for a pair of monochromatic 
photons is required to be around 4 — 8% of the thermal cross section {av) = 3 x 10~^^cm'^/s [2j. 
But, the DM annihilation into monochromatic photon(s) must be loop-suppressed, because 
dark matter is charge-neutral. Thus, it is challenging to find a model for suppressing the tree- 
level annihilation channels into the SM particles while the loop-induced annihilation channels 
are sizable [7-12 . For instance, one needs a large coupling between dark matter and a new 



charged particle running in loops or rely on a resonance effect in the presence of a mediator. 

We consider a singlet Dirac fermion dark matter whose mass can be of weak scale naturally 
due to a new chiral symmetry or gauge symmetry. If there is no direct interaction between 
DM and new charged particles, we need a mediator field connecting dark matter to the SM 
particles through the s-channel. Since the DM annihilation must be s-wave to be a relevant 
process at present, the mediator is either a pseudo-scalar (axion- mediation) (8l[9] or a vector 



[Z' mediation) [10,11 . In this paper, we revisit two representative models for a singlet Dirac 
fermion dark matter with the axion scalar or Z' gauge boson mediators, from the perspective 
of the collider physics. First, we review the mechanisms to generate the Fermi gamma-ray 
line in the considered models: DM annihilation into a pair of monochromatic photons in 
axion mediation and DM annihilation into a single monochromatic photon due to Z7 or /17 
channel in Z' mediation. In axion mediation, a weak-scale axion couples to dark matter and 
also generates U{1) Peccei-Quinn(PQ) anomalies depending on the SM representation and PQ 
charge of a heavy fermion running in triangle loops. In this case, the axion resonance enhances 
the DM annihilation cross section into a photon pair. In Z' mediation, dark matter is charged 
under an extra local U{\) with the Chern-Simons(CS) interactions to the electroweak gauge 
bosons. Thus, depending on the value of the CS coupling, the DM annihilation cross section 
into a single photon can be enhanced with or without a resonance. We recap the consequences 
of the two models in view of the Fermi gamma-ray line, including the relic density condition. 
In particular, we first include a new DM annihilation channel into WW in Z' mediation in the 
presence of an extra CS term for the W-boson. This new channel provides a freedom to take 
the Z^'-boson mass even above twice the DM mass, being consistent with the relic density. 

Searches for a monophoton or monojet with large missing energy at the Large Hadron 
Collider (LHC) have become a powerful constraint on dark matter and large extra dimensions 
|13[[l4| . In the effective operator approach for dark matter interactions with the SM quarks, a 
DM pair can be produced directly from the scattering of an initial quark pair so a monophoton 
radiated by the initial quark is recoiled against the missing energy |15|. In the DM mediator 



models that we consider, however, the effective coupling between dark matter and the SM 



^The preliminary result of the Fermi LAT collaboration has confirmed the 130 GeV gamma-ray line, which 
is shifted to 135 GeV with the reprocessed four-year data though [4]. There is a concern with a bump observed 
in the Earth limb data at the similar energy plls]. 
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quarks is absent after the mediator is integrated out[^ Thus, dark matter cannot be produced 
directly from an initial quark pair, but instead the mediator field is produced in association 
with a monophoton. If the mediator decays into a DM pair, it is possible to impose the 
bound of the monophoton search with large missing energy at the LHC |5 Below the mass 
threshold of a DM pair, we rely on the mediator decay modes into ZZ or Z7 and we expect 
the high efficiency of the photon cut and missing energy if at least one Z — )• z/P is required. 
We consider the collider bounds on the coupling and mass of the mediator field from the 
monophoton searches at the LHC 8 TeV and discuss the prospects at the 14 TeV upgrade. 
In axion mediation, the axion production cross section is suppressed due to a small coupling 
between the axion and the electroweak gauge bosons, so there is no bound at the LHC 8 
TeV. On the other hand, in Z' mediation, if the CS coupling between the Z'-boson and the 
electroweak gauge bosons is sizable, the Z' production cross section is large enough to be 
accessible in the monophoton search already at the LHC 8 TeV. 

The paper is organized as follows. We begin with the description of the dark matter models 
and discuss the constraints on the model parameters from Fermi gamma-ray line and DM relic 
density. Then we present the production cross sections of the mediator fields at the LHC in 
both DM models and set the collider bounds on the models from the monophoton searches 
at the LHC 8 TeV and 14 TeV. Finally, conclusions are drawn. There are two appendices 
showing the detailed computations of the DM annihilation cross sections and decay rates, and 
the mediator production cross sections. 

2 Models for Fermi gamma-ray line 

We consider fermion dark matter models with axion scalar or Z' gauge boson mediators and 
discuss the constraints on the models coming from Fermi gamma-ray line and relic density. 
The stability of fermion dark matter is ensured by the Z2 symmetry, x ~^ which can be 
the remnant of the breaking of PQ symmetry or local U (1). 

2.1 Axion-mediated fermion dark matter 

First, in axion mediation (8||9], after the U{1) PQ symmetry is spontaneously broken, the 
effective action for a Dirac fermion dark matter x ^ pseudo-scalar (axion) a can be 
written in the following, 

^axion = ixi^d^x - m^XX + \{d^.af - - Jf^.F*^'' 

^One can introduce the scalar partner of the axion mediator with Higgs portal coupling in axion mediation [s] 
or the kinetic mixing between Z' and Z bosons in Z' mediation [ll], for which the effective DM couplings 
to the SM quarks appear at tree level. Then, those couplings can be constrained by DM direct detection 
experiments. 

''in the case of a scalar dark matter, the DM annihilation into a photon pair is induced by loops containing a 
charged scalar [7 . In this case, it would be possible to have a direct production of a monophoton in association 
with a scalar DM pair too. 



2 



where ci,C2 are constant parameters given by ci = Ti:{qpQY'^) and C2 = Tr(gpQ/(r)) with 
gPQ(y) being PQ charge (hypercharge) and l{r) being the Dynkin index of representation 
r of the extra lepton under the SU{2)i. When dark matter gets mass by the VEV of the 
scalar partner of the axion, the dark matter couphng is constrained by A^^- = \/2m^/vs- The 
effective axion interactions can be rewritten in terms of physical electroweak gauge bosons, 

-^anomaly — C^y^CL Cf^upcrFl^ Ff^ ~\~ Cz-yd (ifiupaF z 

+czza e^^p^F^^F^" + cwwa e^^paF^ F^ (2) 

where 

c^-y = — ^ — (ciai cos^ 9w + C20i2S:Vc? Ow), (3) 

cz^ = rrrr — (C2a2 - cioi) sin(26'vy), (4) 

czz = z-^ — {c2a2 cos^ Ow + ciaisiT? Ow)-, (5) 

cww = o — • y^) 

We note that the axion mass appears to break the PQ symmetry explicitly, but it can 
come from a PQ invariant higher dimensional operator, such as -^^S"^^^, where a complex 
scalar S = s + ia contains our weak-scale axion as the imaginary part and $ is another 
complex scalar breaking the PQ symmetry spontaneously at high scale by ($) = with 
10^ GeV ^ v$ < 10^^ GeV and consequently containing the invisible axion. In the presence 
of the higher dimensional operator, there is a mixing between the invisible axion a$ and our 
weak-scale axion a^. For the S singlet with PQ charge —2, the higher dimensional operator 
leads to ^ ~ 02/^5 — 2ai/v^ and ^ ~ ai/w^ + 2^302/^^ for ^ Vs, in terms of the mass 
eigenstates 01^2 of invisible and weak-scale axions. Thus, the couplings of a$ to extra heavy 
quarks [16^| or the SM quarks [17| generate the QCD anomalies dominantly by oi, whereas 
the couplings of a<j to extra heavy leptons generate the electroweak anomalies dominantly by 
a2 (9). 

In this model, the DM annihilation cross sections |^ into a photon pair and Z7 are given 
by §[9] 

^ '''^ 27r {Aml-mlY+Tlml' ^' 

^ 4^ (^ml-mlY + Tlml\ Ami) 

where the decay rate of the axion is given by 

Ta = Tain) + Ta{Z^) + T a{Z Z) + Va{WW) + T a{xx) (9) 

with 



r,(77) = ^|c,,|2, (10) 

TT 



*The DM annihilation cross sections for the other channels, WW, ZZ, Z'y, are given in Ref. [s]. 
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Figure 1: Left: Region between nearby black solid lines is the parameter space of rua vs 
Vg satisfying (av)^^/ {av)o = 0.04 — 0.08 where {(tv)q = 3 x lO^^^cm^/s is the thermal cross 
section. Region between nearby blue dashed lines is the parameter space of nia vs Vg producing 
the relic density within WMAP 5a. Right: The same contour plots for nia vs C2/C1. The 
dot-dashed lines are the contours of (o"w)^^/(2(cjf )-Y-y). We have set dark matter mass to 
my = 135 GeV for Fermi gamma-ray line at 135 GeV. 



m" 



L ,0/. Mix 3 



"a 



ra(Z7) = . (11) 



r.,iZZ) = ,12) 



Va{WW) = ^\cww?['^ f , (13) 
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Taixx) = ^mjl-^) . (14) 
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The last decay mode into a dark matter pair is allowed only for > 2m^. When the anomaly 
interactions ([2]) are generated by an extra charged lepton with PQ charge (7pq, electromagnetic 
charge Q, mass rrij and axion coupling A/, the anomaly coefficient for a photon pair is given 
by ci -|- C2 = qpqQ^ in the decoupling limit of the extra charged lepton. We can implement 
the mass dependence of the extra charged lepton with qpQ = 1 and Q = — 1 in eq. ([7]) by 
replacing (ci -|- C2)/vs by ^ r arcsin^(l/-y/T) with r = m'j/m'^ joj. In the limit with r S> 1, 
we can recover the calculation with effective axion interactions. 

In order to explain the observed Fermi gamma ray line by dark matter annihilation, one 
needs the cross section of dark matter annihilating into a pair of monochromatic photons to 
be {av)^^ = (1.27 ± 0.32+°;^^) x lO'^^cm^s"^ for the Einasto profile and {av)yy = (2.27 ± 
0.57lo;5^) X 10"^^cm^s"^ for the NFW profile, that is, Br(xx 77) - 4-8% for thermal dark 
matter [2]. In axion-mediated dark matter models, the DM annihilation into a photon pair 
can give rise to Fermi gamma-ray line at 130(135) GeV for = 130(135) GeV near the axion 
resonance with ~ 2m^, for which the loop factor is compensated for a large annihilation 
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cross section. Then, the DM annihilation into Z7 can produce another gamma-ray hne at 
= 114(120) GeV. The branching fraction of the two-photon channel and the decay width 
of the axion depend on the anomaly coefficients j|8jj . There are microscopic models with heavy 
vector-like leptons for the PQ anomalies where Fermi gamma-ray line and dark matter relic 
density can be explained by axion-mediation at the same time [9]. 

On the left of Fig. [T| we depict the parameter space compatible with the DM annihilation 
cross section for xx ~^ 11 ^'^^ the relic density condition with WMAP 5cr bound, on the vs 
Va plane, for a set of anomaly coefficients with ci = 1, C2 = 3, which is the case for a vector-like 
doublet {Y = —1/2) fermion and a triplet {Y = 0) fermion running in loops |^ On the right of 
Fig. [l| we show the dependence of the Fermi gamma-ray line and relic density conditions on 
the ratio of anomaly coefficients for a fixed singlet VEV, Vg/ci = 55 GeV. In the same figure, 
we also overlay the contour lines of the intensity ratio of one-photon to two-photon lines, 
r = (av) Z'y/i2{(Tv)-y^). We note that when C2/C1 > 1.5, there is a parameter space compatible 
with both Fermi gamma-ray line and relic density, for which the intensity of the one-photon 
line is comparable to or larger than the one of the two-photon line, for instance, r ~ 1 for 
ci = 1, C2 = 3. The presence of two photon lines is an interesting feature of the axion mediated 
DM models. For the collider study in the later section, we have scanned the parameter space 
{vs, ci, C2), in the ranges of 55 GeV < Vs < 200 GeV and < ci,2 < 10, for a given axion mass, 
having in mind the Fermi gamma-ray line and the relic density condition. We note that the 
continuum photons coming from WW, ZZ, Z'j are consistent with the gamma-ray [25j and 



anti-proton |26j constraints as well as the shape of the gamma-ray line spectrum 27 



We also note that when the axion is promoted to a complex scalar in a UV completion, the 
real scalar partner of the axion leads to extra s-channel DM annihilation channels into a pair of 
the SM particles through Higgs portal term In this case, those extra channels contribute to 
the thermal cross section at freezeout and determine the relic density correctly together with 
the axion-mediated channels for general anomaly coefficients, but those channels are p-wave 
suppressed at present. Then, thanks to the extra s-channels with the real scalar partner, 
we are free to choose the anomaly coefficients to be consistent with the Fermi gamma-ray 
line, without taking into account the relic density condition seriously. However, the mixing 
between the real scalar and the Higgs boson can be constrained by the XENONIOO bound 
and the Higgs data and Higgs-like boson search at the LHC (sl. 



2.2 Z'-mediated fermion dark matter 

Second, in Z' mediation, a fermion dark matter has a charge under extra local U{1) and the 
Z' gauge boson has the effective interaction with a single photon only due to Landau- Yang 
theorem. The effective action for fermion dark matter x is 

Cz' = ixl^D^x - m^XX - ^i^^.i^'^" - \inl,Z'^Z'^' + Ccs (15) 



^Extra vector-like leptons can enhance the Higgs-to-diphoton rate by their Yukawa couplings to the Higgs 
doublet [9 24 . The interplay between Fermi gamma-ray line and Higgs diphoton signal in the models with 
extra vector-like leptons has been discussed recently j9] . 
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= 149 GeV 




mz' [GeV] 

Figure 2: Region between solid lines is the parameter space of mz' vs gz' ai satisfying 
{(Tv) z-y/ {c^v)o = 0.08 — 0.16 where {(7v)o = 3 x 10~^^cm^/s is the thermal cross section. We 
have set dark matter mass to my = 149 GeV for Fermi gamma-ray line at 135 GeV. 



where the covariant derivative is -D^x = (dfi — 'i'dz' Z'^/2)x and the effective Z' interactions 
can come from the Chern-Simons(CS) terms as follow^ 

Ccs = a.e'^-^P^Z'^Z^Fj^ + a-^e^^-f- Z'^Z.F'^^ + b^e^^'^P'^ Z'^{W+ d^W- - d.W^W-). (16) 

Here, F'^^ = dpZ'^ — df,Z'p and 01,02,61 are the dimensionless parameters generated by the 
triangle loops with extra heavy fermions. We assume that there is no direct coupling between 
the Z' gauge boson and the light quarks in the SM, apart from a small kinetic mixing with the 
Z-boson. Then, we can avoid the strong dijet bound and the electroweak precision constraint 
for the Z' gauge boson. 

Furthermore, we can write the effective interaction between the Z' gauge boson and the 
Higgs as 

A£c5 = ^6^^'"^Z'a./iFX (17) 



where 03 is a dimensionless parameter and M is the cutoff scale. As in the model in Ref. |llj, 
when top quark couples to the Z' gauge boson, the Higgs coupling 03 can be generated by 
top quark running in the triangle loop. 



The first coupling a\ in eq. (16) gives rise to the DM annihilation into Z7, which can ex- 



plain the Fermi gamma-ray line. Using the results in appendix A, we get the DM annihilation 
cross sections into Z7 and ZZ, respectively, as 



{av)zy 
{(^v)zz 



g^iol cos^ 9y/ 



16m^ 



327r 

4 



(4m^ 



m|,)2 + T\,m\, 
16m^ 
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1 



1 



4m| J 

2 \ 5/2 



~'~ 4m2 



m 



(18) 
(19) 



The first two terms were considered in fermion dark matter models in Ref. 



10 



6 



where the decay rate of the Z' gauge boson is given by 

Vz' = Tz' [Zj] + Tz' (ZZ) + Tz' (WW) + Fz' (xx) (20) 

with 

r.,z.) ^ ||„;c„.^,„.(i-^)'(-L.-L). pi, 

r-<«) = 48!lp<"'-+^'"x)(l-3)' ■ 

We note that when mz' < rri^, the gauge couphng to dark matter opens up the t-channel 
annihilation into Z'Z', with the cross section, 

, V 9%' "^Y / ml, \ 3/2 

'"''>^-^' = & (,.|,-2K)- ('-4) ■ l^'^' 

Furthermore, there are extra DM annihilation channels coming from the other CS cou- 
plings in eq. (16). First, for mz' < 2m^ — raz, the CS coupling 02 contributes the s-channel 
annihilation into Z' Z as 



{av)z'z 



where 



gl,almly'\A,Rz^,Rz) 



1287r 



\{20,Rz,,Rz) _L iRz'+Rz 
4 



+ 4^Rz^ ^(4, Rz',Rz) + 4(8i?z' - Rz) 
{Am^ — m|,)2 + r|,m|, 

(26") 



A(x, y, z) = + + z"^ - 2xy - 2yz - 2zx, (27) 
Rz' = ('^)\ Rz=(^y. (28) 



Second, the CS coupling bi for the W-boson also leads to the s-channel annihilation into a 
W~^W~ pair, with the cross section, 

0-^ w - 2 TTi::? — i:2-^^T^r^^-i:2- 1 - • (29 



647rm^ (4m| — m^,)^ + r|,m|, \ m| y 

In Z'-mediated dark matter models, there is no DM annihilation into a photon pair due 
to Landau- Yang theorem. Instead, the DM annihilation into Zj or /17 leads to a gamma-ray 
line at E-y = Tny,{l — m^/j/(4m^)). Thus, in order to obtain the Fermi gamma-ray line at 
130(135) GeV, we need to take = 144(149) GeV for Z7 or = 155(160) GeV for hj 
where we chose the Higgs mass, nih = 126 GeV. In Fig. [2j we depict the parameter space 
compatible with the DM annihilation cross section for xX ~^ on the mz' vs ^z'^i plane 
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Figure 3: Region between black solid lines: parameter space of mz' vs ai satisfying 
{(yv) z-y / {(yv)^ = 0.08 — 0.16 where {(tv)o = 3 x 10~^^cm^/s is the thermal cross section. 
Region between blue dashed lines: parameter space of m^' vs 6i producing the relic den- 
sity within WMAP 5cr. We have chosen gz' = 0.2 (Left), 0.5 (Right) and set DM mass to 
= 149 GeV for Fermi gamma-ray line at 135 GeV. 



for m-)^ = 149 GeV. In order to get the correct photon flux, we need the DM annihilation 
cross section into a single photon to be twice the fitted value of the DM annihilation cross 
section into a photon pair as discussed in the previous section. 

Now we discuss the constraint of the relic density on the model. First, if < mz' < 
2m^ — mz, there is an extra channel with Z'Z coming from the a2 coupling, which allows 
the possibility to obtain a thermal cross section. In this case, the second coupling 02 in 
eq. (16) is constrained to satisfy the WMAP constraint. Second, if mz' < m-x, the t-channel 
DM exchange leads to the main contribution to dark matter annihilation so the relic density 
condition is independent of a2. As a result, the relic density condition can be satisfied by the 
Z' gauge coupling or the 02 coupling, in the low Z' mass region below — mz- Lastly, for 
mz' > — mz, there are two final states with ZZ and Z^ coming from the ai coupling. 
But, since {crv) z z / Wv) z-y ~ 0.3, the Z7 final state is always dominant. In this case, we 
could not get the correct relic density without extra channels. We note that if kinematically 
allowed, the CS coupling of the W-boson in eq. (16) can open an extra annihilation channel 
with a W^W~ final state. Then, it is possible to explain the relic density with the bi coupling 
only, even for mz' > 2mx — mz- In the following discussion, we focus on the CS couplings 
ai and 61, in particular, in the high Z' mass region. 

We overlay in Fig. [3] the parameter for mz' vs bi satisfying the WMAP bound on the relic 
density within 5a, in comparison to the parameter space of the oi coupling shown in Fig. [2j In 
the low Z' mass region below m^, the contribution of the WW channel to thermal cross section 
can be bounded due to the Z'Z' channel, depending on the value of the Z' gauge coupling. 
On the other hand, for mz' > m^, the thermal cross section is determined dominantly by 
the WW channel, so it only depends on the product of the couplings, gz'bi- We find that 
when the Z'Z' channel is subdominant, the CS coupling for the VF-boson should be larger 
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than the one for the Z-boson by order of magnitude. We note that the continuum photons 
coming from the WW channel and other channels without monophotons are consistent with 



the gamma-ray 25 and anti-proton [26] constraints and can be also compatible with the 
shape of the line spectrum |27|. 

Now we comment on the case that the top quark couples to the Z' gauge boson with the 
same Z' charge as dark matter. In this case, the same Z' gauge coupling is shown to generate 



the CS coupling to the Higgs boson as given in eq. (17). From the 03 coupling in eq. (17), 



we obtain the annihilation cross section for XX ~^ ^7 follows, 

_ gl^gjcos^giy lem^ / \' ...s 

^ ^'^^ 327rM2 (4m| - m|)2 + r|,m|, V 4m2 ) ' ^''^^ 

So, the DM annihilation into /17 can give rise to an extra line at the energy smaller than the 
photon energy at the Fermi line. From eq. ( |18| ), the annihilation cross section for xX ~^ ^7 
is comparable to the one when xX ~^ ^1 for M/a^ ~ mz/ai. In this case, we can get two 
gamma-ray hues at = 116(122) and 130(135) GeV for = 144(149) GeV from /17 and 
Z7, respectively. Furthermore, if kinematically allowed, there are extra decay modes of the 
Z'-boson as follows. 



Ncgl, ,0 o, / 4m? \ 1/2 

^Z'{tt) = — ^^(m|, + 2m2)(l- ^ . 



ASirmz'''"^^ ' ""K ml. 



(33) 



On the other hand, the top quark coupling to Z' can open a new DM annihilation channel 
into a ti final state and allow for Z'-boson to be heavier than 2m^ — mz, in which case the 
ZZ and Z7 annihilation channels only cannot explain the relic density. Suppose that the 
VF-boson CS coupling is small. Then, for mz' ^ 200 GeV and m^ < rrif, the relic density can 
be obtained by the top quark coupling to Z' , when m-^ is slightly above or below mz'/'^^ i.e. 
close to the Z' resonance, or Z' couplings are large, for which slightly smaller than mt is 
preferred. For m^ > m^, the DM annihilation proceeds dominantly into a ti final state, so, for 
moderate couplings, gz' ~ 1, there is a robust prediction for the dark matter mass, whatever 
the value mz' ^ 350 GeV: m^ ~ 150 GeV, which is far away from the Z' resonance 11 . In 
this case, we note that the Z'-boson can also decay into a ti pair with the branching fraction 
larger than the one of the xX decay mode. 



2.3 Summary 

In axion mediation, Fermi gamma-ray line at 130(135) GeV is explained mostly by the DM 
annihilation into a photon pair and it requires m^ = 130(135) GeV and favors 02/ ci > 1 
for the correct branching fraction of the corresponding annihilation cross section. Since we 
rely on the axion resonance for the enhancement of the annihilation cross section, the axion 
mass range is typically constrained to be close to = 260(270) GeV for ci, C2 of order one. 
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Figure 4: The relation between lines in the gamma spectrum and possible signatures at the 
LHC. 

Varying Vg (on the left of Fig. [T]) and/or 02/ ci values (on the right of Fig. [T]), we may relax 
the axion mass up to as large as 500 GeV for Vg ~ 55 GeV and (ci, C2) ~ (1, 10). 

In Z' mediation, the DM annihilation into Z7 or h'y can explain Fermi gamma-ray line 
at 130(135) GeV for = 144(149) GeV or 155(160) GeV. The correct value of the corre- 
sponding annihilation cross section can be obtained for gz'O-i ^ 0.01 or gz'a3 ~ O.OlM/mz, 
far away from the Z' resonance. If there is a VF-boson CS coupling, there is no limit on the 
Z'-boson mass coming from the relic density, because we can vary the M^-boson CS coupling 
to satisfy the relic density condition. In this case, the relic density is obtained by the WW 
annihilation channel with gz'bi ~ 0.1, far away from the Z' resonance. If there is no VF-boson 
coupling, the relic density condition requires mz' < 2m^ — mz in the absence of the top cou- 
pling to Z' and mz' ^ 200 GeV in the presence of the top coupling to Z' . The top quark 
coupling to Z' would open up an extra Z' decay mode into tt when mz' > 2m-^. In order to 
maximize the monophoton events with large missing energy in the next section, henceforth, 
we assume the case with a sizable 1^-boson CS coupling and no or small top gauge coupling 
to Z'. 



3 Mediator production at colliders 

In this section we consider the production of the axion and Z' mediators at the LHC. The 
same diagram which mediates dark matter annihilation into two vector bosons, can be rotated 
to lead to LHC production of the mediator, see Fig. [4j 

Unlike the approaches with the effective interactions between dark matter and the SM 
quarks the dominant production is in association with a vector boson 

pp^V* + X (34) 

where X = a, Z' and V = 7, Z or W^. One could also produce the mediator in vector boson 
fusion (VBF) processes 

pp^X + jj (35) 



where the jets j are rather forward |30|. We will discuss this mode in Sec. 4.3 



The collider searches of a monophoton coming from the initial state radiation have been considered in 
Ref. [T5I. 
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3.1 Axion production 



In axion-mediated dark matter models, the same interactions producing the DM annihilation 
into a photon pair can lead to monophoton events in association with a dark matter pair, 
— )• 7 a, at the LHC. 

When rua > 2m-^ in the half region of the axion resonance, the axion can decay into a dark 
matter pair dominantly for the dark matter coupling A^^ of order one, see eqs. 14 If lUa < 'irn^ 



in the other half region of the axion resonance, the dark matter production is kinematically 
suppressed so the production of a pair of electroweak gauge bosons is dominant. In the case 
where the axion decays into ZZ or Z7, one would obtain missing energy from the decay of 
a Z-boson into a neutrino pair. Here, we assume that the extra vector-like lepton running in 
the loops is heavy enough as < 2m j such that the extra vector- like charged lepton pair is 
not produced from the axion decay. 

For the monophoton production at the LHC, the only relevant parameters are the axion 
decay constant Vs, the anomaly coefficients, 01^2, the axion mass m^ and DM mass m-^. The 
parton- level production cross section for qq ^ is |^ 



aiqq -^07) 



2a 



+ 



+ 



Zy 



4g 



rrir. 



sm{26w) (s - + r|m| 



where v„ 



r|(i 



4|Qg| sin^ 9w) and Uq 



T|, that is. 



1 



1 - ^ sin^ Ow], Vd 



2 V 3 

Similarly, the production cross section for qq 



a{qq aZ) 



2a 



+ 



- 1 + - sin^ 9w 
aZ IS 



-ad 



(36) 



(37) 



.^'"^^^^ sin2(20H/) {s 
4:QqVqCzyCzz s{s - m| 
sm{26w) 



m|)2 + T'^^m? 



is 



At 



z 

1 1 
m„ mi 



(38) 



3.2 Z' production 

In Z'-mediated dark matter models, due to the same CS coupling leading to the DM annihi- 
lation into Z'y, there is a similar process, qq — >• jZ' , where off-shell Z-boson produces a Z' 
gauge boson in association with a monophoton. 

For mz' > 2m^, Z'-boson can decay into a DM pair, producing a large missing energy. In 
order for the branching fraction into XX to be dominant, we need gz' ^ oi, ^ij from eqs. (21 )- 



(24). Thus, from the conditions for Fermi gamma-ray line and relic density obtained in the 
previous section, gz'ai ~ 0.01 and gz'bi ~ 0.1, we get the bounds on the CS couplings, 
ai ^0.1 and 61 < 0.3, for the dominance of the xX decay mode. Below the mass threshold 
of a DM pair, Z'-boson decays into Z7 and/or ZZ,WW. In the case where the Z'-boson 
decays into Z7, we have two photons and missing energy coming from the decay of a Z-boson 



'See Appendix B for details. 
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into a neutrino pair. When the Z'-boson decays into ZZ, we have monophoton events with 
missing energy coming from four neutrinos, but the ZZ decay mode is suppressed by sin^ ^VF; 



see eqs. (24). 



For the monophoton production at the LHC, the relevant parameters are the Z' CS 
couphng a\ and the Z' mass mz' ■ The parton-level production cross section for — >• Z' 7 is 



cT(gg-^Z'7) = , " ^1 + ^ 1- — ) • (39) 



See Appendix B for details. 



4 Setting Limits 

We extract the trivial dependence on couplings by defining a re-scaled cross section 

a = 5 (40) 

where 



c2 



c = I for the axion, and a\ for the Z' , (41) 

Vsjv 

with w = 246 GeV. 

For illustration, we show in Fig. [5] the scaled cross section a in association with a photon. 
The axion cross section is much smaller because in the coefficients Cj one has already pulled out 
the ai/SiT prefactor, plus the axion coupling have one more derivative. This extra-derivative 
makes the off-shell vector boson to have a larger invariant mass [29], resulting in a PDF 
suppression. Nevertheless, if Vs is small and ci,2 relatively large, the model can be accessible, 
see the red-dotted line in Fig. [5} 

In setting limits, we focus on the signature of a photon recoiling against missing energy. 



Both CMS |13| and ATLAS 14 carry this search, where the selection criteria is as follows 



1. Isolated photon with \'q^\ < 1.44(2.37) and > 145(150) GeV (42) 

2. Missing Energy $t > 130(150) GeV (43) 

where the numbers in parenthesis correspond to ATLAS. There are other isolation criteria, 
such as the requirement that the missing energy and hadronic activity, or the photon, should 
not be too close. Veto on a number of high-p^ jets or leptons is also applied, but there is no 
veto on extra photons. For details, look at the experimental notes. 

Although the public notes correspond to the 7 TeV run in 2011, one can estimate the reach 



of the 2012 run at 8 TeV as shown in Ref. 28 . Assuming we combine CMS and ATLAS, we 
expect a limit on the cross section 

{ax BRx Acc) < Ifb (44) 

where BR is the branching fraction of a final state and Acc is the acceptance. 

In the following we are going to first discuss the axion model for the LEP bound and the 



8 TeV run (Sec. 4.1 ) and next the Z' model for the 8 TeV run (Sec. |4.2]). We also discuss the 



expected reach to search for axion and Z' in the 14 TeV upgrade (Sec. 4.3) 
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Figure 5: Parton-level production cross section for qq — )• {a,Z') + 7 at LHC 8TeV. Red- 
dashed line corresponds to the axion case, black-solid is the Z' case. For illustration, we show 
the axion red-dotted line with ci^2 = 10 and Vs = 55 GeV. 



4.1 Bounds on the axion from the 8 TeV run 

There are three ways in which the axion is sensitive to the monophoton search. Above the 
DM threshold, we will consider the process 

PP ~^ ^* /Z* — 7 a 
with a^xx- 

Below the threshold, we will rely on the processes 

PP ~^ 'J* /Z* —7- ZZ'j or Z77 
with at least one Z ^ 1^1? , 



with Z — 7- z^P and/or W ^ ul . 

As shown in Fig. [5j the axion case would be difficult to access with the 8 TeV run, unless 
f s <C f and ci^2 are order one. From the point of view of UV models, Vs is given by the scale 
of the PQ symmetry breaking with a complex scalar S containing the weak-scale axion joj, 
and it is related to the mass of a heavy fermion coupled to S and the mass of the real partner 
of S as 

Vs = V-2^ = ^ (45) 



A 



/ 

with A / being the Yukawa coupling of the heavy fermion to S and As being the singlet quartic 
coupling. In order for the effective theory for dark matter and axion to make sense, we need 
to impose mf,ms > nia ~ 270 GeV, so we get Vs > 54 GeV for As = 47r. One can use this 
bound, and other constraints to obtain a bound on as follows: 
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Figure 6: Upper panel: (Left) Production cross section e+e^ — )• 7a — )• 777 at LEP, 
s/s = 181 — 209 GeV. (Right) Branching fraction of the axion decays below 200 GeV. Lower 
panel: (Left) Cross sections of associated production p,p ^ a, V at the LHC 8 TeV in axion 
mediation. (Right) Production cross sections multiplied by branching fraction for 7+MET+X 
(magenta) at the LHC 8 TeV and Z{-^ l+l') + MET + X (blue) at the LHC 7 TeV and the 
limits for the corresponding processes (red solid and dashed lines, respectively). DM mass 
is set to = 135GeV and parameters {vg, 01,02) are chosen to satisfy the Fermi gamma 
ray ratio and WMAP within 5a from Fig[TJ Input parameters are chosen to give the maxi- 
mum and minimum cross section for each axion mass, denoted as the boundary lines of the 
production cross sections. 



Vs > 54 GeV, ci,2 < 10 , ~ 135 GeV, and {a^^v) ~ lO^^'^cm^s"^ 
^ma< 500 GeV. 

In Fig. § we show in the upper panel the multiphoton production cross section at LEP 
with -y/s = 181 — 209 GeV and the branching fraction of the axion decays in the low axion 
mass region and also present in the lower panel the reaches from monophoton limit at the 
LHC 8 TeV and from mono-Z limit at the LHC 7 TeV f3l||32|. The cross sections for all 
the associated axion productions are computed for parameters {vg, 01,02) in the ranges of 
55 GeV < f s < 200 GeV and < 01^2 ^ 10, satisfying the Fermi gamma ray ratio and WMAP 
within 5a. As seen in the left of the upper panel, the axion production cross section with the 
axion decaying into a photon pair is below 6.2 fb at LEP with ^/s = 181 — 209 GeV, which 
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Figure 7: (Left) Acceptance for pp — )• 7XX ^''^^ other channels involving monophoton + MET 
due to neutrinos as a function of mz'- (Right) Total cross section, a x BR x Acc (in fb), 
normalized to ai = 1. 

is within the experimental uncertainties of the measured cross section around 7 — 9.5 pb at 
LEP [33| . On the other hand, we find that the axion production cross sections are below the 
current sensitivities for both monophoton and mono-Z limits at the LHC. 

Finally, we note that astrophysics also constrains the axion-photon coupling and the small 
axion masses. If axions are produced inside a star and escape, they provide an additional 
cooling channel. Then, the non-standard energy loss prolongs the red giant phase and shortens 
the horizon branch phase. In our case, the bound of star cooling on the axion-photon coupling 
can be satisfied even for the weak-scale Vg , if the axion mass is ^ 300 keV [34] . 

4.2 Bounds on Z' from the 8 TeV run 

In this section we focus on the Z' models to discuss the reach of the 8 TeV run. There are 
two ways in which the Z' is sensitive to the monophoton search. Above the DM threshold, 
we will consider the process 

pp-^ Z* -^-fZ' 
with Z' ^XX- 

Below the threshold, we will rely on the processes 

pp Z* ^ ZZ^ or Zjj 
with at least one Z ^ vv . 

We note that in the presence of the VF-boson CS term, there is an extra process for the Z' 
production in association with the M^-boson and the Z'-boson can decay into a WW pair. 
In this case, the produced VF-boson(s) can also provide missing energy through the leptonic 
decay. 

In the first production mechanism, pj- = $Ti up to radiation. Since ruDM — 150 GeV, 
one would expect large missing energy, and high efficiency on the cuts. We have obtained the 
efficiencies for the ATLAS and CMS cuts for the mass range mz' G [300, 500] GeV, and those 
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Figure 8: (Left) Exclusion on the {mz',ai) parameter space with the overlay of the region 
satisfying {av)z'y/{o'v)o = 0.04 — 0.08 where {av)o = 3 x lO^^^cm^/s is the thermal cross 
section. (Right) Exclusion on the {gz'^^z') parameter space. In both plots, red region is 
excluded. 



are in the 50-65% range. Hence the monophoton searches are ideal to search for Z' in this 
mass range. The blue(upper) line on the left of Fig. [T] shows the acceptance for the — t- 7 xx 
channel as a function of the Z' mass at the LHC 8 TeV. 

In the second production mechanism, one relies on the invisible branching ratio of the Z, 
which is 20 %, and the cut on missing energy and high-p^^ photon may not be that efficient. 
Nevertheless, for events such as pp — )• Z'7 and Z' — )• 7Z the efficiency of the photon cut 
is order one, whereas the missing energy cut efficiency is reduced, as only one of the three 
bosons decay to neutrinos. The total efficiency is in the range 30-50%, plus the reduction due 
to the branching ratio Z — )• viy. Various lines on the left of Fig. [7] shows the acceptances for 
the individual channels giving rise to pp ^ X + ^ + MET as a function of the Z' mass at the 
LHC 8 TeV. 



To obtain the efficiencies we have created new models in Feynrules 19 . We then inter 



faced with MadGraph 20 using the UFO model format f2T]. We incorporated hadronization 



and showering effects using PYTHIA |22] and detector effects with Delphes |23j. Finally, we 
have implemented the ATLAS and CMS cuts within ROOT |18j. In our simulation, jets are 
always anti-kr jets of size R = 0.5. 

From the monophoton searches at the LHC 8 TeV run, on the right of Fig. [TJ we apply the 
acceptance for 7-I-MET channels and show the total cross section (cr x BR x Acc) in fb, with 
re-scaled oi. As a consequence, on the left of Fig. [8| we present the exclusion limit on the 
mz' vs oi plane by assuming that Br(Z' — )• xx) = 1 above the mass threshold of a DM pair 
and overlay the region satisfying {'Jv)z'y/{o'v)o = 0.08 — 0.16 where {(tv)q = 3 x 10~^^cm'^/s 
is the thermal cross section, for two values, gz' =0.1, 0.3. We find that the region away from 
the resonance in Fig. [2] can be excluded depending on the Z' gauge coupling gz>- the excluded 
mass range is mz' < 200 GeV for gz' = 0.3, and mz' < 260 GeV and mz' > 360 GeV for 
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gz' = 0.1. On the right of Fig. [sj we also show the exchided region in the {mz' , gz') parameter 
space being compatible with the Fermi gamma-ray line. 

We remark how the exclusion limit on the cross section changes when the extra decay 
mode of the Z'-boson such as WW exists. As discussed in the previous sections, the relic 
density condition requires extra DM annihilation channels than the one in ZZ, Z7. In the 
low Z' mass region below the threshold, we can use the DM annihilation into a Z' Z' pair 
with the Z' gauge coupling for the correct relic density for mz' > 2?tt,^ — mz and the DM 
annihilation into ZZ' channel with the extra CS coupling 02 for < mz' < 2m^ — mz- 
These extra channels do not change the Z' decay rate. However, in the high Z' mass region 
above the threshold, we need extra channel, for instance, the WW channel. In this case, from 



eqs. (23) and (24), there is an extra contribution to the Z' decay rate, which can be written 

(46) 



as compared to the other dominant channels above the threshold as follows, 
_ Tz'{WW) _ ( biV m%, (1 - 4m2^,/m|,)5/2 



rz'(xx) \9Z' ) m^(m|, + 2m2) (1 - 4m2/m|,) V2 ' 

As a result, when hi ~ gz', the WW decay mode can be comparable to the decay mode into 
a DM pair, reducing the production cross section for pp — )• 7XX above the threshold. But, 
the final states of the WW decay can contain missing energy too and the efficiency of the 
corresponding channel is given by Acc{WW) = 28 — 47% on the left of Fig. [7| as compared to 
Acc{xx) = 50 — 65% for the pp — )• 7XX process. Then, the effective production cross section 
for — )• 7 + MET + X above the threshold is given by 

a.,{pp ^ 7 + MET + X) = aipp ,Z') • ^ccixx) + r Ac^WW) _ 

1 + r 

Therefore, even for r ~ 1, the "effective" acceptance is given by ^{Acc{xx) + Acc{WW)) = 
39 — 56%. Consequently, we still have a high efficiency for the wide range of the parameter 
space where the DM production cross section is reduced. 



4.3 Dedicated searches at 14 TeV 

When the LHC is up again in 2015, it may run at 13.5 TeV, and collect a luminosity per year 
of about 100 fb-^ With the increase of energy, the cross section in associated production 
increases by a factor 0(3 — 10), see Fig.joj Therefore, one could estimate that the sensitivity in 
signal could increase by a factor 0(10) per year, as compared to the 2012 run we discussed in 
the previous section. Keeping in mind that this is an optimistic estimate, it would mean that 
the reach in oi would increase by about factor three, and larger masses will be kinematically 
accessible. 

In the low mass region {mx < 2m;^) a dedicated search would be more suitable. In this 
region, the efficiency to the missing energy cut in the ATLAS and CMS search is relatively 
low, and a bump search in the channel 

pp^V* ^ VX with X ^ 77 , (48) 

would be a more efficient way to look for the mediator X = a,Z' . One could also use 
the vector boson fusion channel, where the background rejection gain by cuts on Ay^-j and 
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Figure 9: Parton-level production cross section for qq — )■ {a,Z') + 7 at LHC 14TeV. Red- 
dashed line corresponds to the axion case, black-sohd is the Z' case. For illustration, we show 
the axion red-dotted line with ci^2 = 10 and Vg = 55 GeV. 

hadronic activity could compensate the low parton luminosities of vector bosons fusing into 
the mediator. 

Finally, let us mention that other channels, such as mono-Z or mono-lepton, could be used 
in association with the monophoton production we have discussed here. For a discussion, see 
Ref. [32]. 



5 Conclusions 

We have considered the bounds of monophoton searches on the models of the Fermi gamma- 
ray lines, focusing on the singlet fermion dark matter with a pseudo-scalar or vector mediator. 
Due to the suppression of the direct couplings between dark matter and the SM light quarks 
in the models, the mediator coupling between dark matter and a pair of electroweak gauge 
bosons gives rise to a dominant production mechanism for the mediator field in association 
with a monophoton at the LHC. Thus, depending on the branching fraction of the invisible 
decay modes of the mediator, we have estimated the limit on the mediator production cross 
section with the LHC 8 TeV run and discussed the implications on the models. 

In axion mediation, where a sizable DM annihilation into a photon pair is obtained by the 
axion resonance, the LHC 8 TeV monophoton searches have not reached the sensitivity yet. 
But, the fact that the axion mass is bounded to about 500 GeV within the resonance band 
could give a helpful guide to the monophoton searches at the LHC 14 TeV run. On the other 
hand, in Z' mediation, the DM annihilation into a single photon can be enhanced without a 
resonance, so a sizeable mediator coupling is allowed. In this case, the monophoton searches 
at the LHC 8 TeV already rules out some of the parameter space in the mediator mass and 
coupling. The dedicated searches for the Z' or axion resonant production at the LHC 14 TeV 
are necessary and complementary to the monophoton searches discussed in this paper. In 
particular, a bump search on the 77 channel is especially promising in the low mass region. 
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Appendix A: DM annihilation cross sections and Z' decay rates 
in Z' mediation 

In this appendix, we present the details of the computation of the DM annihilation cross 
sections and the Z'-boson decay rates in Z' mediation. 



From the Chern-Simons term (16), we compute the dark matter annihilation cross sec- 
tion into Z7. The Feynman rules for the Z'^^ — — Afj{j)2) vertex and x — x ~ -^^ are 
2ai cos Qw^iivpaV2 ^ 5z'7m/2) respectively. Then, the scattering amplitude for x(^2)x(^i) ~^ 
Z{pi)'y{p2) are 

^«-^ = (Ml,^.) (^1.,^^^)^ (A,l) 



where 



M^^^z' = ^gz'vihhMh), (A.2) 

Then, the annihilation cross section is given by 

1 2 

The amplitude squared averaged over the spins of the dark matter is 

X T.'7{P2)ei{p2)f2Pi ■ ' • (A.5) 

The trace part is calculated as 

Tr[(|^2 + rn^)Y'{h " = 4{k^k^' + k^k^') - 4(fci • fcs + m^)^'^'^'. (A.6) 

We note that the polarization sums are 

pf^p^' 



z 
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We also use the formulas for the contracted epsilon tensors as follows, 

eM.pae^V'''''^'5^^'9'^"'5^''' = -6^?', (A.9) 

e^^.p^e^'^'P'^^'g^^'g'"''' = -2(<<5^'-<J^'), (A.IO) 

- + bib^^b^;. (A.ll) 

Using the above formulas, we get 

e^vp.e^'u'p'.'g^^'U''' - ^^^g^^'M = -6pl + ^(pipl - iPi-P2?), (A.12) 

^iii/pa£fj.'u'p'a'{kik2 +k2k^)(g 2~j9 P2P2 



rrir 



'^{{ki ■ k2)pl - {ki ■ P2){k2 ■ P2) 

+ {ki ■ k2){pi ■P2f - (ki ■P2){k2 ■Pi){pi ■P2) 
-{ki ■pi){k2 ■P2){pi ■P2) + {ki ■P2){k2 ■P2)pi 



{ki • k2)plpl + {ki • pi){k2 • pi)p2 



(A.13) 



Therefore, for = 0? we get 



1 



|A^|2 = ^a^COS^^iy-^ 2 ^9 T.2 2 

4 ' (s-m|,)2 + r|,m|, 



is-mlf + 16iki-p2){k2-p2) 



rrir 



ki{pi ■P2) -Pi{ki ■P2)^ {k2{pi ■P2) -Pi{k2 ■P2)^ 



(A.14) 



For non-relativistic dark matter, we get s ~ 4m^, so the squared amplitude becomes 



9z' 



alcos'^ew I, 2 2^2(^'^\ 

— -2— — ^— ^ (4m^ - m^) — 2^ + 1 . 



(A.15) 



Then, including the phase factor by 
get the annihilation cross section as 



M 



[1 — m\/s) \M\'^ and using eq. (A. 4), we 



{av)z-y 



gz'O-i cos 6w 



16m': 



m 



2 \3 



1 



1 



m| ~'~ 4m| 



(A.16) 



327r (4m| - m|)2 + r|,m|, V ^m^, J 

In a similar matter, as given in the text, we can obtain the annihilation cross sections for 
XX ~^ ZZ' , WW from the CS terms in eq. (16) and XX ~^ ^7 from the CS coupling for 
the Higgs in eq. (17). Furthermore, when mz' < m^, there is a t-channel annihilation channel 
for XX ~^ ^' ^' due to the Z' gauge coupling to dark matter. 
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Now we consider the decay rate of the Z' gauge boson. The decay rate is given in terms 
of the decay ampHtude, M-z'^z-^i = ^z'iP^ + P2)-M.'^,^^^, by 

2 



1 



Mz' 



The squared decay amphtude averaged over the spins of the Z' gauge boson is 

cos^ 9w ef^upaef^'v'p'a' ^ {pi + P2)4' (^"1 + ^"2) x 



(A.17) 



E 'ziPlVz iPl) E 'AiP^VA iP2)p'2 pi ■ 



Then, using eqs. (A.9)-(A.ll) and -p^ = 0, we get 

= lal cos^ Ow {ml, - rn^f (\ + -\-) . 



(A.18) 



(A.19) 



Therefore, including the phase factor by 
we get the decay rate as 



Z'^Z-y 



"^alcos'ew (1-5 



l-m|/m|,) and using eq. (|A.17) 
2 \ 3 



1 



1 



,2 + ^2 



(A.20) 



Similarly, as given in the text, we can get the decay rates for Z' — )• ZZ, WW, xx from the CS 
coupling in eq. (16) and for Z' — )• hjjhZjtt from the CS term for the Higgs in eq. (17) and 
top quark coupling to Z' . 



Appendix B: Mediator production cross sections at the LHC 



In this appendix, we present the details on the derivation of the production cross sections 
for the axion and Z' gauge boson in association with a monophoton at the LHC. 

First, in axion mediation, we compute the Drell-Yan production cross section for the 
axion in association with electroweak gauge boson at the LHC: qq' — )■ Vi — )■ aV2. For the 
quark gauge coupling, g/Vi,^/'7^(t'/ - af-f^)f, and the axion coupling, cviVi^^^iupaO. F^^^ F^^ , 
ignoring the quark masses, we get the squared scattering amplitude for q{k2)q' {ki] 
a{pi) V2{p2) as 



V* 



1 



4s 



2„2 



v9fCviV2 



m 



Tr 



£27 



^^fiupa^fi'i/' p'<t' ^ ^ 
V2 



(^V2iP2)^V2(.P2) P3P3 P2P2 



2wjaj7 



where sy is the symmetric factor, which is equal to 1 for Vi 7^ V2 and 2 for Vi 
after calculating the trace and the Lorentz contracted quantities, we obtain 



(B.l) 
V2. Thus, 



2 324 5/C 
N, 



2-2 C,,2 



Pi{ki ■P2){k2 ■P2 



{Pi ■P2)\ {k2 ■Pi){ki ■P2) + (/C2 ■P2){ki -pi) 



P2{k2 ■Pi){ki - Pl) 



(B.2) 
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Consequently, we get the production cross section as follows, 

\P\ 



1 4 9fcl,v,{v} + a}) \p\ 



{s-ml-ml^)(E^E2 + \\p\') 



mv,{El-\\p\' 



where 

\P\ 
El 
E2 

In particular, for 
process is 



(ma + my2)^ 1/2/ ^ma-mv2?\^/^ 



-)'( 



2 
2 

^ 7 



/ rni — m^r \ 



(B.3) 

(B.4) 
(B.5) 
(B.6) 



07, we get sv = 2, so the production cross section for the 



fq (Cl + C2fa^ 



3iVc (167ri;,)2 



(B.7) 



with being the electromagnetic charge of the quark. If the axion decays into a DM pair, 
there are a single photon plus missing energy in the final states. The full cross section for 
a 7 is 



a{qq -^07) 



2a 

3iv: 



sin^(26lu^) {s-m^zf + T^z^^z 



AQqVqCj^czj s{s - m|) 



m 



2x3 



sm{2ew) (s-m|)2 + r|m|_ 
where Vq = T3 (1 — 4:\Qq\ sin^ 0w) and Uq = T§, that is. 



(B.8) 



(B.9) 



Second, in Z' mediation, we compute the Drell-Yan production cross section for Z'{Z) in 
association with a monophoton at the LHC: qq ^ ^ Vi V2. For the quark gauge coupling, 

gqV3^f^qj'^{vq — aq^^)q, and the CS coupling, aie'^'^'^'^ Z'^Z^F^^, ignoring the quark masses, we 
get the squared scattering amplitude for 5(^2)^(^1) ^ ^"3* Vi{pi) 1^2(^2) as 

h7'''kY{v^g + al-2vqaq^'> 



\M\ 



2 _ 



Tr 



1 gqol cos^ 6w 
xe;,.,,e^V'p'a' Yl ^vM)4M) E ^vM)4M)P2Pi- (B-IO) 



V2 
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Thus, after calculating the trace and the Lorentz contracted quantities, we obtain 

pl{ki ■ k2) + 2{ki ■P2){k2 ■P2) 



+ - 



1 



m 



{Pi ■P2)[ik2 ■Pi){ki ■P2) + {k2 ■P2){ki -pi) 



-p\{ki ■P2){k2 ■P2) -pl{k2 ■Pi){ki -pi] 



Consequently, we get the production cross section as follows. 



a{qq Z'7) 



247riVc (s-m|)2 + m|r| 



s\ 1 + 



(B.ll) 



(B.12) 



with Qa = 7i — ^ — . 

When the SM light quarks couple to Z' boson through a kinetic mixing e between Z' 
boson and Z-boson, we get gq = tgz' and Vq = \ and = 0. Then, we also obtain the 
Drell-Yan production cross section for Z in association with a monophoton as 



a[qq ^ Z i) 



1 e^5^,af cos^ Qw — 

m% J \ s 



247riVc (s - m|,)2 + m|,r|, y ,10^ 



2 \ 3 



(B.13) 



In the text, we consider this process to be subdominant, as compared to the Z'7 channel. 
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